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REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM 

OF FROBENIUS 

MATTHIAS BECK AND SHELEMYAHU ZACKS 

Abstract. We study the Frobemus problem: given relatively prime positive integers ai, . . . , ad, 
find the largest value of t (the Frobemus number g{ai, . . . ,ad)) such that '^'l_j^mka.k = t has no 
solution in nonnegative integers mi, . . . ,md- We introduce a method to compute upper bounds 
for g{ai,a2,az), which seem to grow considerably slower than previously known bounds. Our 
computations are based on a formula for the restricted partition function, which involves Dedekind- 
Rademacher sums, and the reciprocity law for these sums. 



1. Introduction 

Given positive integers ai < 02 < • • • < with gcd(ai, . . . , ad) = 1, the linear Diophantine problem 
of Frobenius asks for the largest integer t for which we cannot find nonnegative integers mi, . . . , rrid 
such that 

t = mioi H h nidad . 

We call this largest integer the Frobenius number g{ai, . . . , Od)', its study was initiated in the 19th 
century. One fact which makes this problem attractive is that it can be easily described, for 
example, in terms of coins of denominations oi, . . . , a^; the Frobenius number is the largest amount 
of money which cannot be formed using these coins. For d = 2, it is well known (most probably at 
least since Sylvester [?]) that 

(1) g{ai,a2) = aia2 - ai - a2 . 

For d > 2, all attempts to find explicit formulas have proved elusive. Two excellent survey papers 
on the Frobenius problem are [?] and [?]. 

Our goal is to establish upper bounds for g{ai, . . . , ad)- The literature on such bounds is vast; it 
includes results by Erdos and Graham [?] 

Ol 



(2) 5(01 ) • ■■,ad) < 'iad 
Selmer [?] 

(3) g{ai,...,ad) <2ad-i 
and Vitek [?] 

(4) g{ai,. . . ,ad) < 



d 



Od 

d 



ad 



\{a2-l){ad-2) 



1 . 
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Here oi < 02 < • • • < a^, and [x\ denotes the greatest integer not exceeding x. Davison [?] 
established the lower bound 

(5) 5(01,02,03) > \/3aia2a3 -01-02-03 . 

Experimental data [?] shows that Davison's bound is sharp in the sense that it is very often very 
close to (7(01,02,03). On the other hand, the upper bounds given by ((2), ©, and (|1J) seem to 
be quite large compared to the actual Frobenius numbers. In this paper, we derive a method of 
achieving sharper upper bounds for the Frobenius number. Our results are based on a formula for 
the restricted partition function (Section which involves Dedekind-Rademacher sums, and the 
reciprocity law for these sums (Section The main result is derived in Section 0] computations 
which illustrate our new bounds can be found in Section [SJ 



We focus on the first non-trivial case d = 3; any bound for this case yields a general bound, as 
one can easily see that 5(01, . . . ,0^) < 5(01,02,03) if oi, 02, and 03 are relatively prime. If not 
then we can reduce by one variable at a time: Again by the definition of the Frobenius number, 
5(01, . . . , ad) < 5(01, . . . , ttd-i) if oi, . . . , ttd-i are relatively prime. If not, we can use a formula of 
Brauer and Shockely [?]: If n = gcd (oi, . . . , ad-i) then 

\ /O^l C^d— 1 \ / 

(6) 5(01, ...,ad) =n g{—,..., , o^ + (n - 1) Orf . 

Hence 

/ \ Z*^! 0'd—l\ , \ 

g{ai, ...,ad)<ng{—,..., + (n - 1 o^ . 

\ n n / 



2. The restricted partition function 



We approach the Frobenius problem through the study of the restricted partition function 

P{ai,...,aa}in) =#[{mi,...,md) G'ZiQ: mioi H h m^o^ = n| , 

the number of partitions of n using only oi, . . . , o^ as parts. In view of this function, the Frobenius 
number g{ai, . . . , o^) is the largest integer n such that P{ai,...,ad}(^) = 0- 

In the d = 3 case, we can additionally assume that o = oi, 6 = 02, and = 03 are pairwise relatively 
prime, a simplification due to Johnson's formula [?]: if n = gcd(o, b) then 

g{a, b,c) = n g (-, -, c ) + (n - l)c . 
\n n J 

(This identity is a special case of ©.) 

In the case that o, b, c are pairwise relatively prime. Beck, Diaz, and Robins derived the following 
result for the partition function P{a,fe,c} ["^j Theorem 3]: 

n^n/lll\l/o6c 

+ S-nib, c; a) + S-nic, a; b) + S^n{a, b; c) . 
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Here [?, Equation (14)] 

m=0 ^ ^ ' ' 

where aa~^ = 1 mod c and ((x)) = x—\x\ — 1/2, is a special case of a Dedekind-Rademacher sum; 
we will discuss these sums in the next section. 

To bound the Frobenius number (from above), we need to bound P{a,b,c} (from below), whose only 
nontrivial ingredients are the Dedekind-Rademacher sums. A classical bound for the Dedekind- 
Rademacher sum yielded in [?] the inequality 

g{a, b,c) < - (^^y abc {a + b + c) — a — b — , 

which is of comparable size to the other upper bounds given by @, and (0J). However, we will 
show that one can obtain bounds of smaller magnitude. 



3. Dedekind-Rademacher sums 



The Dedekind-Rademacher sum [?] is defined for a, 6 € Z, x, y G M as 

k + y 



Ria,b;^,y)=T.ii'-^+- 



where 



k=0 



((x)) ifx^Z, 

if X e z. 



Rademacher's sum generalizes the classical Dedekind sum R{a,b; 0,0) [?]. An easy bound for the 
Dedekind-Rademacher sum R{a, b; x, 0) can be obtained through the Cauchy-Schwartz inequality: 
if a and b are relatively prime then 



\R{a, b; x. 



k=0 ^ ^ 



< 



V b 



E 

\k=Q 



(9) 



'6-1 



\k=i ^ 



< 



'k 1 1^^ 



\k=0 



b ^ \ ( b 1 1 
12 ^ 66 J \T2 ~ 4^ 6b 
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(In the third and fourth step we use the periodicity of ((a;))*.) An important property of R{a, b; x, y) 
is Rademacher's reciprocity law [?]: if a and b are relatively prime then 



(10) 
Here 

Q{a,b;x,y) 
where 



R{a, b; x, y) + R{b, a; y, x) = Q{a, b; x, y) . 



1 + J_ + J_ _L 

^ 12 \b ^ ab ^ a) 



if both x,y £ 



i{x)y i{y)y + i i^My) + i;May + bx) + ^^2(2;)) otherwise, 



^2{x) = {x- [x\f -{x-[x\) + 1/6 



denotes the periodic second Bernoulli function. Among other things, this reciprocity law allows us 
to compute R{a, b; x, y) in polynomial time, by means of a Euclidean- type algorithm using the first 
two variables: simply note that we can replace a in R{a, b; x, y) by the least residue of a modulo b. 



To express S in terms of R, we rewrite ((H| as 

c-l 

St{a,b;c) = ^ 

m=C 

Accordingly, 

St{a, b;c)=R i-a-^b, c; -—^,0 ) + 



+ 



a-H 



b-H 



b-H 



if c\t, 
otherwise. 



if c\t, 
otherwise. 



To ease our computations, we bound this as 
(11) 



St{a,b;c)>R^—a ^b,c; 



4. Upper bounds for g{a, b, c) 



To bound St{a, b; c) from below (which yields an upper bound for g(a, b, c)), we use an interplay of 
(jlUj) and @ to obtain a bound for the Dedekind-Rademacher sum corresponding to St, according 
to 1)11(1 . The idea is to reduce the arguments of the Dedekind-Rademacher sum after the application 
of (|10|1 . which means that the bound given by © will be more accurate. To illustrate this, let ci 
be the least nonnegative residue of — a~^5 modulo c. Then 

(12) R[-a-^b,c- —M =i?(ci,c; —M = Q(ci,c; ^, j - ( c, ci; 0, — \ . 



If ci = 1 then the right-hand side can be simplified, as R [c, 1; 0, —^^-^j =0. If ci 7^ 1 then the 

Dedekind-Rademacher sum on the right-hand side of (|12() can be bounded (via ©) sharper then 
the Dedekind-Rademacher sum on the left-hand side. In fact, by a repeated application of (|1U|) . we 
can achieve bounds which are even better. To keep the computations somewhat simple, we apply 
HlOj) once more and illustrate what this process yields in terms of lower bounds for St- Let C2 be 
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the least nonnegative residue of c modulo ci . If C2 = 1 then 



(13) =Q ci,c; ,0 -i? l,ci;0 



= Q(ci,c; — ,0j - Q( l,ci;0 

as R (^ci, 1; =0. If C2 / 1 then (O can be refined as 

i?(-a~^6,c;-^,0^ = Q (^ci, c; 0^ - (^ca, ci; 0, 

(14) =q(^ci,c;-^,0^ -Q (c2,ci;0,-^^ + (^ci, ca; ) . 

The Dedekind-Rademacher sum on the right-hand side can be bounded according to © as 



(15) R{c.,c,;-^-^,0]>J('^^ + ±](^^-U ' 



12 6C2J \12 4 6c2 

We still need to bound Q. 1^2 has a minimum of —1/12 (at x = 1/2) and a maximum of 1/6 (at 
X = 0). These extreme values yield for 

Q ci,c; ,0 



c 



the lower bound 

(16) Q (ci,c; -^,0) >-i + l| + I^-^ = QioUcuc) , 

as well as the upper bound 

a^^t\ Co 1 ci 

< + + — 

c J - 12ci 12c2Ci 12c2 

These inequalities yield the following. 



(17) Q(c2,ci;0, — ) < + — — + -1- = Q^p(c2,ci) . 



Proposition. Suppose a and b are relatively prime to c. Let ci be the least nonnegative residue of 
—a~^b modulo c, and let C2 be the least nonnegative residue of c modulo c\. 

c 1 3 

(i) /fci =1 then St(a,b;c) > 1 . 

\) J I *^ ' ' ^ - 24 6c 4 

(ii) If c, ^ 1 and 02 = I then ^*(«> ^) > + " ^ " ^ " § " ^ 

(iii) // ci / 1 and C2 / 1 then 



S,ia,b;c)>^ + -^ '--^ Jf^ + — ) f^-- + — 

^ ' ' ^ - 12c 12cic 24ci 12ci 12ciC2 12c2 4 V \12 602 J \ 12 4 6c2 
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Proof. (P) Use (O with m = 1 in THl: 



St{a, b;c)>R (^-a-^b, c; 0^ 



Q l,c; 



a-h „\ 1 
2 

a-^t \ 1 



c J 2 
1 1 c 1 
-~4'^6c~24~2' 



Here the last inequahty follows from ()16() . 



(EH) Use (Uni in (HU together with the bounds ^ and ((T7|) : 

a'^t „\ 1 
2 



5t(a, b;c)>R (^-a-\ c; -^,0^ 



Q ci,c; ,0 -Q l,ci;0 



7-1 -VI— I n 

c J \ c J 2 

- 4 12c 12cic 24ci \6ci 12 J 2 



(jml Use (dH) with the bounds given in ^E^, and (fT7|) : 

a-H „\ 1 
2 



5t(a, b]c)>R (^-a-^b, c; 0^ 




a-^t „\ 1 
2 



These lower bounds can be combined with Q and the quadratic formula to give an upper bound 
on the Frobenius number. 

Theorem. Suppose a, b, and c are pairwise relatively prime. Denote the lower bounds for St{b, c; a), 
St{c,a;b), and St{a,b;c) according to the previous proposition by a, (3, and^, respectively. Then 



g{a,b,c) < (a + 6 + c)2 - i (a2 + 62 + - 2abc (q + /? + 7) - i(a + 6 + c). 

One should note that a + /3 + 7 is negative. We can see that the growth behavior of this upper 
bound is dominated by —2abc (a + /3 + 7) under the square root. This means that if we can make 
— (a + /? + 7) somewhat smaller than min(a, b, c) then we get a bound which grows considerably 
less that the bounds given by @, and @. In fact, we can see this difference in example 
computations already when we use the bounds a, (3, 7 as given by our proposition. What is more 
important, however, is the fact that we can easily obtain even better bounds by improving our 
proposition through additional applications of Rademacher's reciprocity law (|1U() . We illustrate 
this with the following algorithm, whose result is a bound on St{a, b; c), which can be used in the 
above theorem (instead of the bounds coming from the proposition). 
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Algorithm. Input: a,b,c (pairwise relatively prime) and N (number of iterations). Output: lower 
bound S for St{a, 6; c). 

c_l := - a~{-l} b modulo c (least nonnegative residue) 



c_2 := c modulo c_l (least nonnegative residue) 
S_l := S + Q_low (c_l,c) 
S_2 := S_l - Q_up (c_2,c_l) 

IF c_l = 1 THEN S := S_l 
Else S := S_2 

IF c_l = 1 OR c_2 = 1 OR n = N THEN BREAK 
c := c_2 

c_l := c_l modulo c_2 (least nonnegative residue) 
n := n + 1 



IF c_l > 1 AND c_2 > 2 THEN S := S - sqrt((c_2/12 + 1/(6 c_2) - 1/4) (c_2/12 + 1/(6 c_2))) 
S := S - 1/2 

The algorithm repeats the steps described in the proposition N times, at each step bounding Q 
coming from Rademacher reciprocity according to (|16|) and (|17j) . It stops prematurely if one of the 
variables is 1, in which case the remaining Dedekind-Rademacher sum is zero. 



In the present section we illustrate the newly proposed upper bound for g{a, b, c) numerically. In 
order to compare the results also with the lower bound given by Davison we present here the 
values 

/(a, b, c) = g{a, b, c) + a + b + c . 
For these Frobenius numbers, Davison's lower bound is 

f{a,b,c) > \/3z , 



where z = V abc. In [?] we presented together with David Einstein an algorithm for the exact 
computation of f{a,b,c). Einstein computed 20000 "admissible" (see [?]) values of f{a,b,c) for 
relatively prime arguments chosen at random from the set {3, . . . ,750}. In [?] we arrived at the 



empirical conjecture that f[a,b,c) < V abc . The objectives of our current presentations are: 



S := 
n := 1 
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